Velocity renormalization and anomalous quasiparticle dispersion in extrinsic graphene 
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Using many-body diagrammatic perturbation theory we consider carrier density- and substrate- 
dependent many-body renormalization of doped or gated graphene induced by Coulombic electron- 
electron interaction effects. We quantitatively calculate the many-body spectral function, the renor- 
malized quasiparticle energy dispersion, and the renormalized graphene velocity using the leading- 
order self-energy in the dynamically screened Coulomb interaction within the ring diagram approx- 
imation. We predict experimentally detectable many-body signatures, which are enhanced as the 
carrier density and the substrate dielectric constant are reduced, finding an intriguing instability 
in the graphene excitation spectrum at low wave vectors where interaction completely destroys all 
particle-like features of the noninteracting linear dispersion. We also make experimentally rele- 
vant quantitative predictions about the carrier density and wave-vector dependence of graphene 
velocity renormalization induced by electron-electron interaction. We compare on-shell and off-shell 
self-energy approximations within the ring diagram approximation, finding a substantial quantita- 
tive difference between their predicted velocity renormalization corrections in spite of the generally 
weak-coupling nature of interaction in graphene. 



I. INTRODUCTION 

Graphene is a unique electronic material which consists 
simply of a single two-dimensional (2D) atomic mem- 
brane of carbon atoms in a honeycomb lattice. The sem- 
inal discovery^ that graphene can be gated (i.e., "doped" 
with a finite free carrier density) to produce a variable 
carrier density (n) electron or hole 2D system simply 
by tuning an external gate voltage has given rise to a 
truly amazing level of intellectual activity with more that 
10,000 graphene publications appearing in the scientific 
literature since 2005 covering many disparate fields in- 
cluding physics, chemistry, materials science, electrical 
engineering, and chemical engineering. From the funda- 
mental perspective of quantum condensed matter physics 
graphene is an intriguing system^ in many ways: it is a 
true 2D stable crystal with unique elastic and mechani- 
cal properties, and gated graphene is an exotic 2D gapless 
and chiral electron system with linear energy dispersion, 
i.e., the electron- hole bands in graphene have Dirac-Weyl 
massless linear dispersion instead of the usual parabolic 
energy dispersion one sees in most solid state materi- 
als. Graphene is thus both a semimetal and a gapless 
semiconductor, which can be easily doped by an external 
gate voltage to induce a variable carrier density n (either 
electrons or holes) with a special charge neutrality point 
(n = 0), often called the Dirac point, separating electron 
bands from the hole bands. 

The purpose of the current work is to theoretically 
consider electron-electron interaction effects in gated 
graphene (i.e. n 7^ 0), sometimes also called "extrin- 
sic" graphene^ to contrast with very special situation 
of "intrinsic" graphene (n = 0) where the chemical po- 
tential is precisely located at the Dirac point. Since the 
Dirac point itself (n = precisely) is by definition a set 
of measure zero (particularly because of the gaplessness 



of graphene bands) , most experiments can only probe ex- 
trinsic graphene (with n ^ 0) with any conclusion about 
the Dirac point being obtained indirectly through an ex- 
trapolation to zero density. The physics of graphene 
Dirac point is interesting in its own right as it is well- 
established^"— that intrinsic graphene (n = 0) is not a 
Fermi liquid. The presence of disorder-induced electron- 
hole puddles^ makes it difficult to study the physics of 
pristine Dirac point physics. Interaction physics at the 
Dirac point is outside the scope of the current work. 

Since our interest is extrinsic gated graphene at T = 0, 
we consider our starting noninteracting system to have 
the Fermi energy (Ep, or equivalently the chemical po- 
tential fi) in the conduction band (with no loss of gen- 
erality) given by Ep — wo(47rn/g)^/^, choosing h = I 
throughout, where g = A is the graphene ground state 
degeneracy arising from spin (2) and valley (2) degenera- 
cies. The Fermi wave vector kp = (tttt,)^/^ defines the 2D 
momentum upto which the conduction band with linear 
noninteracting band dispersion Eo{k) = vok, with k = |k| 
being the 2D wave vector, is filled. The linear band 
dispersion leads to a simple linear E density of states, 
which then provides kp, Ep, etc. quoted abovo^. Here 
fo ~ 10* cm/s is the noninteracting graphene velocity ob- 
tained, for example, from band structure calculations^. 
Electron-electron interaction would affect the graphene 
noninteracting dispersion i?o(fc) = vok, perhaps modify- 
ing the dispersion itself, i.e., changing the constant vq, 
which is independent of both the carrier density and the 
substrate material supporting the graphene layer in the 
noninteracting single-particle approximation, to a renor- 
malized wave vector-dependent nonlinear graphene veloc- 
ity v*(k) which depends not only on k, but also on the 
carrier density n and the substrate material (through the 
background dielectric constant k). 

The main goal of our work is to make precise (and 



experimentally verifiable) quantitative predictions about 
the interacting quasiparticle energy dispersion and the 
renormalized interacting velocity v* in gated extrinsic 
graphene which, in principle, may depend on wave vector, 
carrier density, and the substrate material. Quite surpris- 
ingly, this issue, in spite of its considerable importance, 
has not been studied in details in the existing theoretical 
literature as earlier work on the subject has concentrated 
mostly on interaction effects in intrinsic graphene^"— and 
the associated instabilities of the Dirac point or on the 
numerical calculations of the interacting spectral func- 
tion in a narrow range of momentum and frequency^i^. 
In addition, we investigate the renormalized quasiparticle 
energy dispersion for different interaction strength (i.e. 
Ts) values to investigate instability at low energy regime 
(i.e., near the Dirac point). It is expected that the quasi- 
particle feature may be modified due to the many body 
effects (plasmonic feature) as seen in the ARPESi^rJ^. 
However, unexpectedly we find that the quasiparticle fea- 
tures disappear as the interaction strength increases. In 
earlier worki^ some of the features presented here are dis- 
cussed for small Vg values {vg < 0.5). However, the insta- 
bility due to many body effect is significantly enhanced 
for large r^ and for low densities, which is not discussed 
elsewhere. This is a very interesting prediction and can 
be observed in samples with high r^ value (for example, 
suspended graphene). In fact, we find that a careful ex- 
perimental study of even the standard graphene on sub- 
strates should be able to discern our predicted signature 
of a dispersion instability in the spectral function. 

We emphasize that there have been several theoreti- 
cal investigations of interaction effects in graphene^riSii^ 
including some by us^i^i^ as well as a few careful experi- 
mental studiesiiiSrJ^ comparing with the theoretical pre- 
dictions. Our work is based on exist ing^r^ii^iS and well- 
established theoretical techniques (namely, Feynman- 
Dyson diagrammatic perturbation theory using the dy- 
namically screened Coulomb interaction as the effective 
coupling constant - a technique that goes back to the 
1960s and has been used extensively and successfully to 
calculate many-body interaction effects on single-particle 
properties in metals, doped semiconductors, and 2D elec- 
tron systems as well as in band structure calculations 
under the name of ^GW approximation'). In spite of the 
existing theoretical literature dealing with the calcula- 
tion of electron-electron interaction-induced many-body 
renormalization effects in doped graphene, the results 
presented in the current work are completely new and 
are directly relevant to experimental measurements. Al- 
though based on the standard GW-typc theoretical ap- 
proximation our numerical results for the coupling con- 
stant dependent graphene Fermi velocity as well as the 
wave vector dependent graphene velocity are unavail- 
able in the literature. These quantities can be measured 
directly by STM/STS or ARPES measurements to be 
compared with our theoretical results. Similarly, our 
results for the anomalous quasiparticle dispersion and 
the associated collapse of the quasiparticle picture for 
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FIG. 1. Self-energy contribution [(a)-(d) plus all higher-order 
ring diagrams] included in our theory with the wiggly and 
the non- wiggly lines representing the bare interaction {V) and 
Green's function (Go) respectively. Figures (e) and (f) rep- 
resent two second-order contributions ignored in the theory, 
which may approximately cancel each other leading to the 
on-shell approximation being quantitatively "better" than the 
off-shell Dyson's solution. 



graphene at low carrier densities are new although the 
graphene spectral function itself has earlier been calcu- 
lated in the literatureSr"— . Finally, our comparison be- 
tween on-shell and off-shell self-energy approximations 
for graphene many-body renormalization has not earlier 
been investigated in the graphene literature where all ear- 
lier work specifically considered the off-shell solution of 
the Dyson equation. 

The rest of the paper is organized as follows. In Sec. 
II we discuss the general theory of the self-energy of 
graphene. In Sec. HI we present our calculated renor- 
malized Fermi velocities within on-shell and off-shell ap- 
proximations. In Sec. IV we show the detail results of 
the quasiparticle spectral function and the instability of 
the graphene dispersion, and finally we conclude in Sec. 
V. 



II. SELF-ENERGY 

The central theoretical quantity we calculate is the 
graphene dynamical self-energy function S(A:, E), which, 
in general, depends both on momentum (wave vector) 
and energy (frequency) independently. An exact the- 
oretical evaluation of the self-energy for electrons in- 
teracting via the long-range Coulomb interaction (V) 
is, of course, impossible for arbitrary carrier density 
and interaction strength, but an excellent approxima- 
tion scheme, which goes by the various nomenclatures 
such as the GPF-approximation or the dynamical RPA, 
is well-established for parabolic metallic 2D and 3D elec- 
tron systems. This approximation involves calculating 
the self-energy in an infinite order diagrammatic pertur- 
bation theory including all diagrams containing (to all 
orders) the bare Coulomb interaction, the bare electron 
Green's function, and the bare polarizability (i.e., the 
bubble or the ring diagrams formed by the closed loop 
of an electron and a hole propagator) . Equivalently, this 
infinite series of ring diagrams (we show in Fig. [T] the 
self-energy diagrams of the theory) corresponds to the 
leading-order perturbation expansion in the dynamically 
screened Coulomb interaction (W) where the screening 



is done by the electronic dynamical dielectric function 
(e) obtained in the random phase approximation (RPA) 
involving the infinite series of ring diagrams^^. 

The validity of considering RPA type diagrams is based 
on the infrared divergence of the ring diagrams. All non- 
RPA diagrams are negligible compared with the same 
order RPA diagram. The divergence can be controlled 
by the RPA re-summation as shown in Fig. [T] The jus- 
tification for why non-RPA diag rams in graphene can be 
neglected is discussed in Rcf. |lj. We also emphasize that 
there is no known controlled technique for systematically 
calculating quasiparticle properties in interacting elec- 
tron systems other than the GW approximation based 
on RPA. Within the leading-order dynamical RPA the 
self-energy Yj{k^E) of graphene is given byi^ii^ 



j:s{k,iUJn} = ^^X! X! Go,s'{k + q,iLLln +Wn) 

xW{q,lv^)Fss'{k,k + q), (1) 



s' q,ii^n 



where /3 = 1/fcsT is the inverse temperature, W is 
the screened Coulomb interaction, Go{k,iuJn) = [iw„ — 
£'o(fc) + ^J]^^ is the bare Green's function (where a;„, z/„, 
and /i are Matsubara fermion, boson frequencies and 
chemical potential, respectively)^. Due to the chiral 
property of graphene the self-energy in Eq. ^ has an ad- 
ditional term (i.e., -F^s') which does not exist in the non- 
chiral systems. The F^gi arises from the overlap of \sk) 
and \s'k') and is given by Fj^' (fc, fc') = ^ (1 + ss' cos 6*^.^'), 
where s,s' = ±1 are band indices of graphene and Okk' is 
the angle between fc, k' . The screened Coulomb interac- 
tion W{q,Wn) = F(g)/e(g,if„), where V{q) = 2'Ke'^ / nq 
is the bare Coulomb potential (k =background dielectric 
constant), and e(g,w„) is the 2D dynamical dielectric 
function. In the RPA e{q,Wn) = 1 ^ ^(9)110(9, Ji^n), 
where Ho ~ GqGo is the noninteracting polarizabilityi^. 
The expression for S, W, Hq involve multidimensional 
complex integrals over internal momentum and energy. 

Since there has been substantial earlier work in the lit- 
erature on the theory of the electron self-energy within 
the dynamical RPA GW approximation in systems as 
different as 3D electron gas^^, 2D electron gasi^, ID 
electron gasi^, monolayer graphene^i^ i^^i^^ , and bilayer 
graphene^S, we refrain from giving more theoretical de- 
tails on our calculations, concentrating instead on the re- 
sults and their experimental implications for observable 
many-body effects in gated graphene. 

We note that {i) the calculation of the self-energy S 
leads immediately to the interacting Green's function G, 



G- 



Gn 



El which contains all information 



about the single-particle properties of graphene, and that 
(ii) once the bare polarizability Ho, which was earlier 
obtainedi^, is known, the calculation of self-energy [Eq.[T] 
would lead to G since Go and V are known by definition. 
We note that we are utilizing the leading-order "GqW^" 
rather than the self-consistent "GW" approximation in 
our theory. We believe the GqW approximation to be 
more meaningful from a diagrammatic perturbativc sense 



since it is a true leading-order expansion in the effective 
RPA screened dynamical interaction W whereas the iter- 
ative GW approximation mixes in unwanted higher-order 
terms in W which are inconsistent with the leading-order 
expansion in W. 

We emphasize that the graphene self-energy calcula- 
tion is highly nontrivial and extremely subtle because 
one must include both intraband and interband contribu- 
tions in the theory. Remembering that the noninteract- 
ing Green's function Gq = [E — i?o(fc)]~^ has its poles at 
the noninteracting energy Eq = v^k, we get the following 
equation for the poles of the interacting Green's func- 
tion E = Eq + 'E{k,E), which defines a general integral 
equation (Dyson equation) for the interacting graphene 
dispersion: 



E{k) =Eo{k)+Re'E[k,E{k)]. 



(2) 



We note that Yi{k,E) is, in general, complex and the 
quasiparticle energy (damping) for the interacting system 
is given by the real (imaginary) part of the self-energy. 
We solve the integral equation defined by Eq. ^ numer- 
ically iteratively to obtain the interacting quasiparticle 
dispersion E{k) — v\{k)k, and refer to this iterative so- 
lution as the "off-shell" approximation which alludes to 
the full solution of the Dyson equation. An alternative 
approximation for the self-energy, the so-called on-shell 
approximation, goes back to the early days of diagram- 
matic many-body theory^!, and has, in particular, been 
emphasized by Rice^^. In the on-shell approximation, 
one simply takes the first iterative solution of Eq. ([2]) to 
get 



E{k) ^Eo{k)+Rc^k,Eo{k)], 



(3) 



which docs not involve solving the Dyson integral equa- 
tion. This would define an on-shell graphene veloc- 
ity given by E{k) = V2{k)k. Note that the imaginary 
part of on-shell self energy near Fermi surface behaves 
as Im[I][fc, £'o(fc)] ex [£'o(fc) — A^]^ ln[i5o(fe) — /x], where 
/i = Eo{kp)i^ Thus the on-shell imaginary part of the 
self-energy certainly vanishes as fc — > fc^?. In fact, the in- 
elastic scattering rate, an extensively used physical quan- 
tity of much experimental interest, is precisely the on- 
shell approximation to the Im[E(fc, E)]. We will compare 
the on-shell and the off-shell graphene self-energy solu- 
tion within our dynamical RPA scheme because the issue 
of which of these two is a better approximation to the 
Coulomb self-energy problem has been much discussed 
in the literature over the last 50 years^lr— . These past 
discussions of course focused only on regular nonchiral 
parabolic electron systems with only intraband contri- 
butions and our current work is on chiral, linearly dis- 
persing, gapless graphene carriers where both intra- and 
inter-band contributions to the self-energy must be in- 
cluded in the theory (as we do in our current work). 
We note that there is no deep principle that could make 
the " correct" theoretical choice between the on-shell and 
the off-shell approximation within the GqW self-energy 



scheme since both schemes obey the Ward identity and 
the Dyson equation up to the order in the interaction the 
theory is meant to be vahd. (i.e., leading order in the 
dynamically screened interaction W). The only effective 
choice between these approximations is an operational 
one based on the detailed quantitative comparison with 
the experimental data as has already been discussed in 
the past^"^. 

It may be worthwhile for us to make some remarks on 
the on-shell versus off-shell many-body self-energy ap- 
proximation in the context of graphene physics, and why 
this may be a question of substantial importance. This 
issue arises because the dynamical Coulomb self-energy 
can only be calculated approximately in the leading-order 
GW-RPA approximation. If an exact self-energy function 
E(fc, E) is available, then obviously the full Dyson equa- 
tion solution would be necessary to obtain the quasipar- 
ticle energy dispersion E(k), making the ofF-shell veloc- 
ity the only meaningful quantity. It has, however, been 
pointed out^"— that the leading-order screened inter- 
action expansion (i.e., the ring diagram approximation) 
for the calculated RPA self-energy necessitates using the 
leading-order iterative solution of E(fc,i?), i.e., using the 
on-shell approximation E(fc, EQ{k)), instead of the full so- 
lution of the integral Dyson equation so as not to mix var- 
ious orders of approximation in the theory. In particular, 
it has been arguecl22i that the on-shell approximation of 
the self-energy is, in fact, quantitatively better than the 
off-shell approximation within the GW-RPA scheme of 
the Coulomb self-energy evaluation for metallic electrons. 
The subject has been controversial^^— with no definite 
consensus in the literature in spite of the almost 50-year 
history of the topic. Our careful quantitative comparison 
of on-shell and off-shell approximations, when compared 
with the available experimental data in graphene should 
lead to a resolution of this old and important theoret- 
ical controversy. Our results tend to indicate that the 
off-shell approximation works better for graphene as de- 
scribed below. 

We emphasize, so that there is no misunderstanding 
of this point, that if the self-energy is evaluated exactly 
then one must always solve the Dyson equation, and as 
such the off-shell theory defined by Eq. ([2]) is the only op- 
tion for the quasiparticle dispersion in any exact theory. 
The question arises, however, whether solving the Dyson 
equation is always preferable for obtaining the quasipar- 
ticle energy dispersion E(k) even when the self-energy 
has been calculated in an approximate theory. In partic- 
ular, if the self-energy is calculated in the leading-order 
expansion in the dynamically screened Coulomb interac- 
tion, as it is in the current RPA-ring diagram GW theory 
and in all the existing theories in the literature, then it 
is unclear which approximation, on-shell or off-shell, is a 
better quantitative approximation. 

Iterating Eq. ([2]) formally order by order, i.e., solving 
the implicit integral equation for the off-shell self-energy 
defined by Eq. ([2|) through the standard iterative tech- 



nique, we get 

Ei{k)=Eo{k) + i:{k,Eo{k)) 

E2{k) = Eo{k) + I](fc, Eoik) + E(fc,£o(fc))) 

E^ik) = Eoik) + I](fc, Eo{k) + E(fc,^o(fc) + E(fc,-Bo(fc)))) 

Ei{k) = ... 

We note that the leading-order iterative solution above is 
the same as the on-shell approximation defined in Eq. ^ . 
We now note that the GW approximation used in the cur- 
rent (and many other) work involves (see Fig. [1]) calcu- 
lating the self-energy formally as a leading-order expres- 
sion in the screened Coulomb interaction T, '^ W, and as 
such all iterations of the self-energy beyond the on-shell 
approximation must necessarily involve terms which are 
formally higher order in W , which would be inconsistent 
with the leading order GW calculation. This is in fact 
the argument behind the claim that the on-shell approxi- 
mation may be a better approximation than the off-shell 
one within the leading-order GW theory. Within this 
leading-order theory there is no obvious way to claim that 
the full solution of the Dyson equation (i.e. the off-shell 
approximation) is necessary better than the leading-order 
solution of the Dyson equation (i.e. the on-shell approx- 
imation) since the former necessarily includes only an 
uncontrolled subset of higher-order terms in the screened 
interaction whereas the latter is manifestly leading-order 
in the screened interaction everywhere. We therefore pro- 
vide results for both on-shell and off-shell approximations 
for future comparison with experiments. 



III. RENORMALIZED FERMI VELOCITY 

For a given self-energy we can calculate the Fermi 
velocity by differentiating the self-energy with respect 
to the wave vector at k — kp, that is, v*(k) = 
dE[k) / dk\k=kF ■ The renormalized Fermi velocities v*p 
corresponding to the appropriate self-energies of Eqs. (2) 
and (3), respectively, are given by 
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and v*p = v\{k = kp)^ Vp^ = V2{k = kp). In our cal- 
culation of the on-shell velocity at the Fermi surface ( 
k — kp and i?o(A:) ~ Ep), we take the limit of both k 
and w = Eo{k) going to the Fermi surface simultaneously. 
Along this path the imaginary part of the self energy be- 
comes identically zero at the Fermi surface. However, 
taking k and Egi^k) to the Fermi surface independently 
gives rise to the different results and would be incorrect. 
We note that, while Vp (both for on-shell and off-shell ap- 
proximations) depends on the coupling constant (i.e., the 
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FIG. 2. Calculated renormalized Fermi velocity (a) within 
the off-shell approximation, Vp_^/vo, and (b) within the on- 
shell approximation, vp,^/vo, as a function of Vs for several 
different densities, n = 10^", 10", 10^^, and 10^'' cm"^ (from 
top to bottom). The dashed lines represent the asymptotic re- 
sults for small Vs given by Eq. ((6]). In (c) and (d) the renormal- 
ized Fermi velocityies (thick lines) within the off-shell approx- 
imation and within the on-shell approximation are shown, re- 
spectively, as a function of carrier density for several different 
Vs, Vs = 0.5 (dod-dashed lines), Vs = 0.8 (dashed lines), and 
Vs — 2.2 (solid lines). The thin lines represent the asymptotic 
results for small r^ given by Eq. ((6|. 



background dielectric constant k) and the carrier density 
(n), the general renormalizicd graphene velocity v* (fc) de- 
pends also on the momentum k, implying that the inter- 
acting quasiparticle graphene energy dispersion E{k) is 
no longer simply linear in k (as the noninteracting energy 
Eq = vok is) even for a given substrate and fixed den- 
sity. Thus, we have v*{k) = v*{k\n,K); Vp = Vp(ri,K), 
and depending on the approximation scheme (off-shell 
or on-shell) we have two distinct renormalized velocities 
vl 2- A main goal of our work is to obtain the {k, n; k) 
dependence of the graphene quasiparticle velocity. 

Before presenting our full numerical results (Figs. [2] 
and [3]) which is one of our two main results in this work, 
we give the leading-order analytic result for the Fermi 
velocity Vp = v*{k = kp) for gated graphene 
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where r^ = e^ /{vqk) is the so-called (background di- 
electric constant dependent) graphene fine-structure con- 
stant which defines the dimensionless electron-electron 
interaction strength or the Coulomb coupling constant 
for the problem and fcc ~ 1/a (where a is the graphene 
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FIG. 3. The calculated v*{k) as a function of k/kp are 
shown for (a) r^ = 0.4, (b) r^ = 0.88, and rs = 2.2 and a fixed 
carrier density n = 10^^ cm^^. The black dashed (red solid) 
line indicates the renormalized Fermi velocity within off-shell 
(on-shell) approximation. In (c) the velocity within off-shell 
is not defined for k/kp < 0.3 (i.e., there is no solution in 
Dyson's equation correpsonding to the quasiparticle.) In (d), 
(e), and (f) we show v*{k) for carrier density n — 10^" cm^^ 
and ra = 0.4, Va = 0.88, and Vs = 2.2, respectively. Since 
there are no solutions in Dyson's equation correponding to 
the quasiparticle the renormalized off-shell velocity at small 
k is not shown in (d), (e), and (f). 



lattice size) is the ultraviolet cut-off for the problem. 
Throughout this paper kc = 1/(2. 46A) is used. Here 
K = (ks + l)/2 where Kg is the lattice dielectric constant 
of the substrate material. The analytic formula given in 
Eq. ^ has earlier been derived in the literature^, and 
we provide it here for the sake of completeness and com- 
parison. 

It is important to emphasize that Eq. 1^ applies to 
both on-shell and off-shell approximations as they have 
exactly the same leading-order expansions in Ts, and 
Eq. ^ provides an exact expression in the Ts — ^ limit 
since the ring diagrams are the most divergent diagrams 
in this leading-order (in 7'^) limit. A curious feature of 
Eq. ^, which is completely unique to graphene with its 
linear massless and gapless energy dispersion, is that the 
velocity renormalization depends on two distinct dimen- 
sionless coupling parameters, r^ [ex {vok)~^] and kc/kp 
[ex {y/na)~^] arising respectively from the Coulomb in- 
teraction strength (r^) and the ultraviolet lattice cut- 
off (kc/kp). In ordinary parabolic electron systems, the 
Coulomb coupling strength r^ itself depends on the car- 
rier density (e.g., in 2D parabolic electron systems with 
a mass to, r^ ^ m/^Jn)^ and no ultraviolet cut-off is 
necessary. Thus, graphene is a very special and inter- 
esting Coulomb system with the interaction strength r^ 
itself being density independent, but renormalization ef- 
fects being logarithmically density dependent through 
the ultraviolet divergence. Such logarithmic dependence 
on the ultraviolet cut-off, while being common in rela- 
tivistic field theories, is rare in solid state physics where 
the infrared Coulomb divergence is typically the main 



problem. We also mention that the second term (the 
square bracket) in Eq. ^ and the last term [containing 
ln(fcc/fci?)] arise respectively from intraband and inter- 
band transitions. Depending on the carrier density and 
the background dielectric constant, one or the other con- 
tribution could dominate, but in general both must be 
kept in the theory for extrinsic graphene. We note that 
the ultraviolet logarithmically divergent term in Eq. [S] 
goes as \n{n/nc) where ric is the cut-off density defin- 
ing kc , and as such the ultraviolet renormalization of the 
graphene coupling constant can be directly studied by 
tuning the carrier density n with its effect being larger 
for larger Vs. 

In Fig. [2] we show our full numerical calculations for the 
renormalized graphene velocity in the off-shell [Fig. [2Ja)] 
and on-shell [Fig. [2jb)] approximations as a function of 
Tg for different carrier densities (n) comparing the results 
with the analytic formula given by Eq. ([6]). In Figs.[2][(c) 
and (d) the renormalized velocities as a function of den- 
sity n for different r^ values are shown. Several features 
stand out in Fig. [2] (i) the velocity renormalization is 
a strong function of both r^ and n in both approxima- 
tion schemes; (ii) the renormalization is substantially 
stronger in the on-shell than in the off-shell approxima- 
tion; (in) the leading-order analytic formula [Eq. ([6|)] 
applies only for r^ < 0.4 — 0.8 depending on the density 
and the approximation scheme. 

In Fig. [3] we show our numerically calculated v*{k) for 
three different r^ values as a function of momentum for 
the fixed carrier density n = 10^^ cm~^ for both approxi- 
mation schemes. There are several interesting features of 
Fig. [3l (i) both approximations imply weak nonlinearity 
in the renormalized graphene quasiparticle dispersion ex- 
cept for k < OAkp] (ii) the many-body renormalization 
is substantially stronger for the on-shell approximation; 
(iii) at low momentum, k < O.Sfc^, the off-shell (on- 
shell) renormalization indicates nonlinear negative (posi- 
tive) velocity corrections in graphene; (iv) as r^ increases 
or at low densities the renormalized off-shell velocity is 
strongly nonlinear; (w) for large r^ values or at low den- 
sities the renormalized off-shell velocity is not defined at 
small fc, that is, there is no solution in Dyson's equation 
[see Eq. ([2|)] corresponding to the quasiparticle energy. 
Thus the quasiparticle feature at low momentum disap- 
pears at high Vs values or at low densities. We emphasize 
that some aspects of this physics were already apparent in 
Ref. [iOj , but the effect should be much more pronounced 
at higher (lower) r^ (n) values. 

One feature of Figs. [2] and |3] requires particular empha- 
sis since it is quite non-obvious: On-shell and off-shell 
approximations start deviating from each other substan- 
tially already in the rather weak-coupling situation of 
r^ > 0.5. Thus, it is clearly important to know by com- 
paring with experiments which one is the better approx- 
imation. 




FIG. 4. The renormalized quasiparticle energy dispersions 
as a function of wave vector for the electron density n = 10^^ 
cm^^ and for various Ts, (a) r^ = 0.4, (b) rs = 0.88, and (c) 
Vs = 2.2. In (d)(e)(f) the renormalized quasiparticle energy 
dispersions are shown as a function of wave vector for a fixed 
r-s = 0.88 and for different electron densities (d) n = 10^^ 



% (e) n = 10'' cm-^ and (f) n = 10'" cm- 



Green 



dashed (black solid) lines represent the non-interacting (on- 
shell) energy dispersion, and symbols are the solutions of the 
Dyson's integral equation. The dots (triangles) represent the 
dispersion of off-shell quasiparticle mode (plasmaron modes). 
The squares indicate the dispersion of the overdamped mode, 
which has insignificant spectral weight. 



IV. QUASIPARTICLE DISPERSION AND 
SPECTRAL FUNCTION 



Although Figs. [21 and 13] present one of our main quan- 
titative theoretical results showing the density, coupling 
constant, and momentum dependence of the renormal- 
ized graphene quasiparticle velocity, we now provide 
some calculated results for the many-body quasiparti- 
cle dispersion and spectral function since the quasipar- 
ticle velocity is derived from the quasiparticle energy 
dispersion [c.f. E{k) = v*{k)k] which is contained in 
the interacting spectral function of the system. Even 
though wc have used the well known many body diagram- 
matic perturbation theory [GqW approximation) wc find 




A 


h 


=0.4 

A 


J 


(d) 

k=0.9k,. 


- J 


ij 


w 


k=0.5k^. 






J 


M 


kpO.lk, 




-="'(f) 

k=0.95k,, 
k=0.9ku. ^ — /j 


^^ 


I 




k=O.StE^ ■J 






'V 

^ k 


k:-0.7k,r 


^^ 


-0.5k 


k=0.3k,, 


M.lk,, 



-5 -4 -3 -2-10 I 

E/Ep 



FIG. 5. Calculated quasiparticle spectral function A{k, E) 
as a function of energy for different wave vectors and three 
different r^ values. The lines for different wave vectors are 
shifted upward for clarity. In (a), (b), and (c) the density 
n = lQ^^cm~^ is used, and the density IO^^ctti"^ is used in 
(d), (e), and (f). 



very interesting and experimentally relevant results in 
the graphene energy dispersion which have not been dis- 
cussed in the existing literature. These results are quite 
unanticipated. 

In Fig. m we show our calculated quasiparticle en- 
ergy dispersion for different densities and r^ values. In 
(a), (b), and (c) results for a fixed density n — 10^^ 
cm~^ and for different r^ values {vg — 0.4, 0.88, and 
2.2 respectively) are shown, and in (d), (e), and (f) re- 
sults for a fixed r, = 0.88 and for different densities 



(n = 10 



13 



10 and n = 10 cm , respectively) are 



shown in both on-shell and off-shell approximations com- 
paring them with the noninteracting energy dispersion 
Ea{k) = vok. The Fermi velocity in each case is given 
by the k = kp point. We note that r^ = 0.4, 0.88, 2.2 
correspond to graphene on SiC (or BN), Si02, vacuum, 
respectively. The important thing to note is that the 
on-shell approximation indicates very strong quasiparti- 
cle dispersion renormalization, particularly for suspended 
graphene (r^ = 2.2), which may be unphysical, thus per- 
haps ruling out the validity of the on-shell approximation. 
Consistent with the results presented in Fig. [21 the off- 
shell approximation predicts rather small many-body dis- 



persion (and velocity) renormalization, but it does make 
a very dramatic prediction about the complete failure 
of the quasiparticle picture for small momentum in sus- 
pended graphene [rg = 2.2). In particular, the off-shell 

cm~^) 



10^ 



quasiparticle mode for r^ = 2.2 (and n 
completely disappears for k < OAkp (we have checked 
the same to be true for other carrier densities also) , leav- 
ing only a weak plasmaron mode which shows up for all 
Tg values at low wave vectors. Thus a clear prediction 
of our theory is that the quasiparticle dispersion in sus- 
pended graphene (r^ = 2.2) will look dramatically dif- 
ferent for k < 0.5kp with the complete disappearance of 
any particle-like feature in the graphene dispersion. As 
shown in Fig. |4l(f) we find that the similar feature (i.e. 
the complete disappearance of a particle-like feature) de- 
velops as the density decreases. For rg = 0.88 and high 



enough densities [n > 10 



11^ 



^ see Fig. \^h) and (d)] 



the quasiparticle-like dispersion is well defined even at 



(7 = 0, but for n < W 



the quasiparticle feature 



starts to disappear near q ~ 0. Thus, there are two dis- 
tinct ways of approaching 'strong-coupling' in graphene 
either by going to smaller kp (i.e. lower density) keeping 
Tg fixed or by keeping kp (or density) fixed and increas- 
ing Tg. However, these two strong-coupling approaches 
give non-equivalent results as shown in Fig. U) Note that 
in ordinary electron 2D gas only changing rg is allowed to 
approach the strong coupling limit whereas in graphene 
one has one more adjustable parameter to approach the 
strong coupling limit, i.e., by changing Vg through back- 
ground dielectric constant of substrate and by changing 
kp through the back gate voltage. Since some aspects of 
this interesting physics were already observed in Ref. [l0[ 
we believe that experiments at higher (lower) rg (n) val- 
ues should lead to a striking confirmation of our theoret- 
ical predictions. 

This prediction can be directly tested in ARPES^StJ^ 
and STS^"— measurements. Earlier theoretical work^ 
on the graphene spectral function completely missed this 
bizarre collapse of the quasiparticle picture for suspended 
graphene since they^i^ focused on graphene on substrates. 
We note that Vg = 2.2 is a relatively weak-coupling 
regime compared with 3D metals (r^ ~ 3 — 6) and 2D 
semiconductors {vg '^ 5 — 30) where no such dispersion 
instability in the interacting spectral function has ever 
been predicted. 

To further emphasize the low wave vector interaction- 
driven instability in the graphene quasiparticle spec- 
trum discussed above, we show in Fig. [5] our calcu- 
lated interacting graphene spectral function given by 
A{k,E) = —2IinG{k,E), which directly provides the 
spectral strength of E{k) features in the graphene spec- 
trum with the noninteracting result being precisely 
Ao{k,E) = 2TrS[E - Eo{k)]. The calculated spectral 
function in Fig. [S] clearly shows that for Vg = 2.2, the 
single particle-like feature disappears completely from 
the graphene spectral function for k < 0.5kp (0.7kp) 
for ?7 = 10^^ (10^°) cm~^, and the spectral weight basi- 
cally gets distributed broadly over the incoherent back- 
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FIG. 6. Density plot of the quasiparticle spectral function 
A{k, E) as a function of energy E/ Ep and wave vectors k/kp 
for three different r^ values, (a) and (d) rs = 0.4, (b) and (e) 
Ts = 0.88, (c) and (f) r^ = 2.2. In (a), (b), and (c) the density 
n = lO^^cm^^ is used, and the density lO^^cm"'^ is used in 
(d), (e), and (f). 



ground. For r^ = 0.4 (graphene on SiC or BN), how- 
ever, some sharp spectral features persist down to low 
wave vectors, whereas graphene on Si02 {ts = 0.88) is 
intermediate in nature with sharp particle-like spectral 
features being present in the spectrum at low wave vec- 
tors and at high density, but disappearing at low den- 
sity. Our theory, which has direct observable conse- 
quences for graphene experiments, thus predicts that the 
whole linear dispersing noninteracting energy band pic- 
ture in graphene should completely collapse in suspended 
graphene for k < O.Bkp in the 10^*^ — 10^^ cm~^ density 
range where all the spectral weight should disappear into 
a broad incoherent background due to electron-electron 
interaction. We emphasize that earlier graphene spectral 
function calculations^ completely missed this dramatic 
behavior arising from the linear chiral nature of graphene. 

To clearly see the disappearance (instability) of the 
graphene energy dispersion at small wave vectors (near 
Dirac point), in Fig. IHlwe show the calculated spectral 
function, A{k,E), as a color plot in k/kp and E/Ep 
space for two different densities (a),(b),(c) n — IQ^^cm? 
and (d),(e),(f) n = lO^^cm^, respectively. In this fig- 
ure the narrow bright area indicates the highest spectral 



weight. For rs = 0.4 [(a) and (d)] the single particle- 
like linear dispersion has the most spectral weight and 
well defined plasmaron features appear below the sin- 
gle particle-like peak. However for Vg = 2.2 [(c) and (f)] 
there are no well defined peaks at small wave vectors and 
the spectral weight is distributed broadly over the inco- 
herent background. The broadening of spectral weight at 
small wave vectors is more pronounced at low density (f). 
The disappearance (instability) develops gradually as the 
wave vector decreases, and the anomalous dispersion is 
stronger for lower density and higher r^. We suggest 
that spectral function measurements via ARPES or STS 
be carried out on suspended graphene to directly observe 
our predicted instability. 



V. CONCLUSION 

To conclude, we have provided the detailed quanti- 
tative results for nonlinear interacting graphene disper- 
sion as a function of density, coupling constant, and mo- 
mentum, predicting in the process a spectacular collapse 
of the quasiparticle picture in suspended graphene for 
k < O.bkp, which can be directly tested in experiments. 
Our results indicate that many-body rcnormalization of 
graphene velocity is possible up to a factor of two as a 
function of either carrier density or wave vector if a large 
range of density and wave vector can be explored in fu- 
ture experiments. Our work is a generalization of earlier 
works^riS in the literature to large r^ values and lower 
carrier densities where the anomalous dispersion proper- 
ties of graphene should lead to a striking collapse of the 
quasiparticle picture as shown in our results. 

Even though we have used the well known many 
body diagrammatic perturbation theory (GqW approxi- 
mation) our main findings (i.e., Fermi velocity rcnormal- 
ization in both on-shell and off-shell approximations and 
an instability in the graphene energy dispersion) are to- 
tally unexpected. Especially, the disappearance of the 
quasiparticle features at small wave vectors as the inter- 
action strength increases is totally unanticipated and has 
not earlier been predicted or discussed anywhere in the 
existing literature to the best of our knowledge. Since 
the current available experimentsiSrJ^i^Sr^ do not indi- 
cate any signature of this instability we believe that this 
work is important and of interest because the results pre- 
sented in this paper are relevant for future experiments. 
The current available experiments have mostly been per- 
formed at small r^ values for graphene on substrates. 
So, to see our prediction of the instability it is required 
to do experiments with samples having small electron 
density and high r^ value such as suspended graphene. 
Since there are two distinct ways of approaching 'strong- 
coupling' in graphene either by keeping density n fixed 
and increasing rg or by going to smaller n keeping r^ 
fixed, the similar many body fcaturcn may appear as the 
density decreases for a fixed r^ value. Thus, in stead of 
increasing the rs value we can measure these many body 



effects by reducing density in graphene on Si02 or SiC 
substates. 
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